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Abstract 



We present a phenomenological analysis of data on both inclusive and semi-inclusive 
spin asymmetries. We examine the impact of the semi-inclusive results presented by 
SMC on the determination of polarized parton distributions performing global fits with 
different sets of observables. We discuss the flavour dependence of the polarized sea inside 
a nucleon. 
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In recent years a number of theoretical attempts J3j 0? S B S S 0> B @' EU ^° 
determine the polarized quark parton distributions in the nucleon have been performed. 
The deep inelastic polarized structure functions gi(x, Q 2 ) or the asymmetries measured in 
inclusive processes [|ll], [12], [13], [14|, [l^, [16], |17|, [18|, [19], |2(| ^T] are used in phenomenological 



analyses. Such an analysis of the first moment of the structure function = /q 1 g±(x)dx 



pointed out that quarks carry little of the spin of the nucleon |13|, [22], |23|. The reasonable 
suggestion is that the sea quarks and/or gluons are polarized. However, inclusive deep 
inelastic scattering does not provide sufficient information about the flavour separation of 
the polarized sea. Hence different combination of the polarized parton distributions have 
to be measured in order to get more information about flavour structure of the polarized 
sea. 

The measurement of the semi-inclusive spin asymmetries for positively and negatively 
charged hadrons from deep inelastic scattering of polarized muons on polarized protons 
and deuterons provides additional data on required observables. Presently available semi- 
inclusive results [p4|, p5| , p6| can be used to determine polarized valence and non-strange 
sea quark distributions, independently from totally inclusive data. The aim of the paper 
is to combine two kinds of existing data, inclusive and semi-inclusive, to extract polarized 
parton distributions. 

Measurement of the inclusive deep inelastic lepton nucleon scattering gives information 
about the spin asymmetry J27j: 



which in leading order QCD parton model is given by: 

A"( X Q 2 ) ~ E ^ A9(x ' g2) (2) 
l( ' g) - Z q e 2 qq (x,Q 2 ) ' {2) 

where e q is the charge of the g-flavoured quark, q and Aq denotes unpolarized and po- 
larized quark distributions respectively, where q = u,d, s, u, d, s. This is the consequence 
of the fact that in LO QCD 2g x = J2 q e lAq, ~F% = X^e 2 - 1 ? an d of the Callan-Gross 
relation F 2 = 1xF\. First attempt to improve such a model is to use R N (x,Q 2 ) = 
[F 2 N — 2xF-^)/2xF^ ^ 0, which is the ratio of the absorption cross-sections for vir- 
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tual longitudinal and transverse photons (R = o"l/o"t) f28[- In calculations, we use the 



parametrization of R described in Ref.|[29||, which is analogous to the one given in Ref.p 



but fitted to the enlarged set of data on R with new experimental values |3T], £32], [33|, |34 . 
This correction leads to the expression: 

for proton and neutron target (N = p,n). The parton distributions are those of the 
proton whereas for neutron are obtained by the isospin interchange u «-> d. For deuteron 
target case one has to multiply above expression by additional factor 1 — 3/2pD where pu 
is a probability of D-state in deuteron wave function (jpr, = 0.05 ± 0.01) |p5f . 



Analogously, for the semi-inclusive asymmetries, the expression in the same order can 
be written as: 

A Nh ( n 2^ ^ Iz dz 9i (x,z,Q ) 
1 {X ' q) z~ f z dzFf\x,z,Q*y W 
where h denotes the hadron detected in the final state and the variable z is given by 

Eh/E^(l — x) with energies given in 7*p CM frame. The region Z is determined by 

kinematical cuts in measurement of the asymmetries. Summing over positively charged 

hadrons, i.e. 7r + , K + and p, and negatively charged (n~ , K~, p) respectively, we get: 

A N+ ^(x O 2 ) ~ ^M^^(x,Q 2 )Djm (1 r n (x o2)) (5) 
1 {X ' Q) - ejgfo Q 2 )D*m (1 + R { ^ Q )} " (5) 

Here D^(Q 2 ) = f 02 dzD^(z, Q 2 ) and D^(z, Q 2 ) is the fragmentation function which repre- 
sents the probability that a struck quark with a flavour q fragments into a hadron h. To 
reduce the number of independent fragmentation functions one can use charge invariance 
and isospin rotation symmetry as well as assumption for the unfavoured and favoured 



fragmentation |23|, |25j. Further assumption concerning the strange quark fragmentation 
function (e.g. Df + + Df~ = 2D^ + ) reduces the number of independent fragmentation 
functions to 6. Finally the set of different weights in eq. @ is: 

Y,D h u = J2 D u = DZ + +D« + + Di , J2 D u = J2 D u = DZ+D«- + Dl , 

h+ h- h- h+ 

h+ h~ h- h+ 

E D s+H D h s = 2(Df + Df + HQ . (6) 



The presence of different J2h Dg in eq. (Q) enables to examine combination of the polarized 
parton distributions different than in the inclusive case. 



To compare theoretical predictions of eq. (Q) and eq. (|5|) with experimental results we 
have to construct or choose the set of unpolarized and polarized quark parton distribution 
functions. These functions are combinations of the elementary ones, i.e. density of quarks 
with spin parallel to the nucleon spin q + (x, Q 2 ) and density of quarks with spin anti- 
parallel to the nucleon spin q~(x,Q 2 ). In details: q(x,Q 2 ) = q + (x,Q 2 ) + q~(x,Q 2 ) and 
Aq(x, Q 2 ) = q + (x, Q 2 ) — q~(x, Q 2 ). Our assumption is that distributions q + and q~ have 
the same functional behaviour, so there is the only difference in the numerical coefficients 
||. It is not necessarily true for q and Aq because the appropriate coefficients in q + and 
q~ could be equal (or have the same absolute value but opposite sign) and in this case 
equivalent coefficients in q (Aq) vanish. The idea is to use formulas for the unpolarized 
quark parton distributions as an input, then to extract from them formulas for q + and 
q~ distributions just by splitting the numerical constants. 

Previously this idea was explored in Ref.||[], where the latest version of the MRS 
|36| parametrization was used. To test the dependence of final results on the input 
parametrization we have chosen the latest version of GRV parametrization for unpo- 
larized parton distributions |57]].This parametrization gives for the valence quarks at 
Q2 = 4GeV 2 : 

u v (x) = 3.221x~ a436 (l-x) 3 - 726 (l- 0.689x°- 2 + 2.254x + 1.261xi), 

d v (x) = 0.507s- a624 (l-x) 4 - 476 (l + 1.615x a553 + 3.651a; + 1.3a;f), (7) 

whereas for the sea anti-quarks: 

s(x) = 0.0034a;- 1 (l -a;) 6 - 166 (l - 2.392 v^ + T.Og^e 2 - 592 ^^-)- 1 - 15 , 

x 

S(x) = x-\l-xf ™ ' ~~ 



0.00285e 2 003 V ll ^+ (8) 

x ai58 (0.738 - 0.981x + 1.063x 2 )(ln -f mi 

x 

-0.596/1 ™N8.621/i , n A A i ™0.876 



5(x) = 0.107x" a59b (l - x) 8 b2i (l + 0.441x u '^ b + 18.721s 



x It is not the first paper on this subject (see references therein) but the data on semi-inclusive 
asymmetries was included into the fit for the first time. 



where S(x) = d(x) + u(x) is the non-strange singlet contribution to the sea and 5(x) = 
d(x) — u(x) is the isovector non-strange part of the quark sea. For the unpolarized gluon 
distribution we get: 

G(x) = x-\l - x) 5 - 566 [o.0527e 2 - 141 v /nr i + x a731 (5.11 - 1.204a; - 1.91Lr 2 )(ln -)-° A72 . 

(9) 

Generally the unpolarized parton distribution for the valence quarks and the isovector 
non-strange part of the sea can be written in the form q(x) = x aq (l — x) l3q W(x). In other 
cases there is a similar part (1 — x)^, which describes asymptotic behaviour for x tending 
to 1 but terms responsible for behaviour for x tending to are more complicated. 

Now we split the above distribution functions between q + and q~ in order to get 
polarized parton distributions as Aq(x) = q + (x) — q~(x). The asymptotic behaviour 
for x — > 1 (i.e. the value of f3 q ) is the same for all distributions like in the unpolarized 
case and for x — > (the value of a q ) it remains unchanged for valence quarks and the 
isovector part of quark sea. We must be more careful in treating the strange sea and the 
isoscalar part. Assuming that the polarized structure function g x have to be integrable 
one has to split appropriate numerical constants in such a manner that non-integrable 
terms of unpolarized parton distributions disappear in polarized parton distributions (i.e. 
one has to split these coefficients equally between q + and q~). This procedure, of course, 
changes asymptotic behaviour but functions remain integrable despite singular behaviour 
at x — > 0. Our expressions for Aq(x) are: 

Au v (x) = x- 0A36 (l-x) 3 - 726 (A u + B u x - 2 + C u x + D u x 3 z), 
Ad v (x) = x-°- 624 (l-x) 4A76 (A d + B d x°- 553 + C d x + D d xi), 

As(x) = x-^il-xf-^iAs + B^e 2 - 592 ^^-)- 1 - 15 , 

x 

AS(x) = x- 0M2 (l-x) 6 - 356 (A s + B s x + C s x 2 )(\n-) 0m7 , (10) 

x 

A5(x) = x-°- 596 (l-x) 8S21 (As + B 5 x - 876 ), 
where we have introduced 15 new parameters. 

The polarized parton distribution functions must satisfy positivity constraint, 

Aq(x,Q 2 ) <q(x,Q 2 ), (11) 



which leads to several constraints on coefficients in each distribution. Furthermore we 
fix the normalization of the non-singlet distributions using the experimental value of the 
axial charge: 

A<? 8 = 3F — D (12) 

where F and D are the antisymmetric and symmetric SU(3) coupling constants of hyperon 
beta decays |38], j39|. Aq denotes the first moment, i.e. the total polarization of each quark 



(or combination of quarks), which is defined as: 

Aq = [ dxAq(x). (13) 
Jo 

The SU (3) flavour non-singlet combinations are defined by: 

Aq 8 = (Au + Au) + (Ad + Ad) - 2 (As + As), 

Ag 3 = (Au + Au) - (Ad + Ad). (14) 

Assuming that the sea contribution for quarks and anti-quarks are equal, first moments 
of above non-singlet combinations become: 

Ag 8 = Au v + Ad v + 2 AS - 4As, 

Ag 3 = Au v - Ad v - 2A8. (15) 
As we do not fix the first moment Ag 3 we are able to test the Bjorken sum rule ED 



Aq 3 = F + D. (16) 

We do not put A8(x, Q 2 ) = (we distinguish Au and Ad), thus we are able to test 
SU(2) isospin breaking effects. Other first moments that can be calculated using the ob- 
tained integrated quark polarizations are: 

AS = Au v + Ad v + 2 As + 2 AS , 

r? = I Au v + ±-Ad v + ±As + ^-AS - \a5 , (17) 
9 18 9 18 o 

r? = — Au v + - Ad v + -As + — AS + -AS , 
1 18 9 9 18 6 



The remaining 16 coefficients of eqs. fllPf) are determined by fitting the available data 
on the inclusive spin asymmetries for proton, neutron and deuteron targets and on the 
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semi-inclusive spin asymmetries for the proton and deuteron target. The fit is performed 
assuming that the spin asymmetries do not depend on Q 2 . Although the latter assumption 
is not consistent with theoretical predictions (Q 2 -evolution of the numerator of eqs. ([l],|]) 
differs from Q 2 -evolution of the denominator due to different polarized and unpolarized 
splitting functions), it is consistent with experimental observation [11-19]. 

The results for the parameters in eq. (|H]) derived from the fit to data on inclusive and 
semi-inclusive spin asymmetries are presented below: 

A u = 0.175 , B u = 0.301 , C u = 2.010 , D u = 6.752 , 

A d = -0.381 , B d = 0.083 , C d = 0.046 , D d = -2.944 , 

A s = -0.00052 , B s = -0.007 , (18) 

A s = 0.026 , B s = 0.574 , C s = -1.863 , 

A 5 = 0.002 , B s = -0.289 , 




0.0 
-0.2 
-0.4 



-0.6 
0.01 



j^n exp. 

j^ti total 

j^ri incl. 

j^n semi— incl. 




Figure 1: The inclusive spin asymme- 
tries A± total obtained from the total fit 
to the inclusive and semi-inclusive data, 
compared to all existing inclusive data 
points. mcL and A? semi ~ incL are pre- 
dictions which come out from the fit 
to data on inclusive and semi-inclusive 
spin asymmetries, correspondingly. 



6 



For the fit we get x 2 — 147 for 157 d.o.f. (we use 123 data points from totally inclusive 
experiments and 48 data points from semi-inclusive experiment), hence x 2 l^d.o.f. — 0.94. 
Results for the inclusive asymmetries for the proton, neutron and deuteron target are 
compared to the experimental data in Fig. [I]. The comparison of the fitted semi-inclusive 
spin asymmetries for production of positively and negatively charged hadrons from the 
proton and deuteron target to experimental points is given in Fig. 0. Polarized quark 
distributions are presented in Fig. [3[ 




-0.2 I ' ' — ' ' ' ' -0.4 I ' ' — ' ' ' ' 

0.001 0.010 0.100 1.000 0.001 0.010 0.100 1.000 



Figure 2: The semi-inclusive spin asymmetries obtained from the fit to the data on semi- 
inclusive spin asymmetries, compared to recent results presented by SMC [2(\] . total 
denotes semi-inclusive asymmetries obtained from the total fit to inclusive and semi- 
inclusive data. Predictions for semi-inclusive asymmetries calculated using distributions 
which come out from the fit to inclusive data are also presented (A^ ). Note the last 
data point for A\ + , which gives the largest contribution to x 2 , even for the semi- inclusive 
fit. 
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Figure 3: The distributions derived from 
the total fit and from fits to the data on 
inclusive and semi-inclusive spin asym- 
metries separately. 
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The values of the first moments of parton distributions are as follows: 
Au v = 0.60 ±0.01 , Au = 0.08 ±0.02 , 

Ad v = -0.56 ±0.01 , Ad = 0.07 ±0.02 , (19) 
As = -0.042 ± 0.004 . 

From these numbers one can evaluate the values of the first moments of the structure 
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functions and other combinations of the polarized quark parton distributions: 

Am = 0.68 ±0.02 , Ad = -0.49 ± 0.02 , 

rf = 0.142 ±0.002 , 17 = -0.057 ±0.005 , , . 

AS = 0.26 ±0.01 ,_ Ag 3 = 1.19 ±0.07, ^ 
Aq sea = 2(Au + Ad + As) = 0.22 ± 0.04 . 

In our model quark contribution to the spin of the proton is dominated by the sea polar- 
ization. The contribution of each valence quark is almost the same but has the opposite 
sign, hence valence quarks carry little of the spin of the proton. However all distribu- 
tions are fitted to the data points from the measured region, i.e. for x > 0.003, and 
contributions beyond the measured region (J ' 003 dx f(x)) are questionable. Moreover the 
low x behaviour of the polarized quark parton distributions is determined by the un- 
polarized ones, therefore it is not consistent with the Regge theory prediction. Though 
application of Regge theory is incompatible with the pQCD HI], E2| it is interesting to 
compare predictions of the Regge type behaviour for low x to the predictions of our 
model, especially as some of the considered quantities change rapidly for x < 0.003. 
For example Jo 003 dx Ad v (x) = —0.113 while / 1 003 dx Ad v {x) = —0.451. It is the conse- 
quence of the power-like behaviour of the type x~ 0,624 . More dramatical change can be 
observed for the non-strange sea quark distribution. The contribution of the low x region 
is J° m3 dxAS(x) = 0.072 while for x > 0.003 we get fl 003 dxAS(x) = 0.079. This is 
due to the fact, that AS(x) is more singular for x tending to than any valence quark 
distribution. Multiplication of the term x" 842 by the term (In ^) 037 enlarges the value of 
integral over the range from x = to x = 0.003, but not significantly. The strange quark 
distribution behaviour for low x is the most complex one. The term x~ 0,5 is suppressed 
with the (ln^) -1 ' 15 term but the multiplying term e 2 - 592 V ln * increases the contribution 
to the integral of the low x region very fast. Finally an integration over x below 0.003 
gives —0.012 to the total polarization of the quark s which is —0.043. 

Now we can compare our results to the more stable Regge theory prediction. The 
quantities integrated over region from x = 0.003 to x = 1 (integration over the region 
covered by the experimental data plus extrapolation for higher x) extrapolated to x — 0, 
postulating the Regge type behaviour for all quark parton distributions of a type x -0 5 , 
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give: 

Au v = 0.61 , Ad v = -0.54 , Au = 0.06 , Ad = 0.05 , As = -0.04 , , , 

r? = 0.132 , T™ = -0.063 , AS = 0.20 , Ag 3 = 1.17 , Aq sea = 0.13 . ^ ' 

Our results on first moments of the proton and neutron structure functions are in 



agreement with experimental results given in Ref . [yj| . Other estimations in Ref . |Tj, [15], 
2"T| are slightly smaller but our results are still consistent within two standard deviations. 
For the first moment of the deuteron structure function gf = | (gf + g" ) (1 — 1.5po) we get 
= 0.039 ±0.004 (or Tf = 0.032 if the Regge behaviour for low x is assumed) what is in 



excellent agreement with results in Ref. [16], yj. For the purely non-singlet combination of 
the structure functions (gf — g™), which in our model is Ti~ n = |Ag 3 , we obtain Ti~ n = 
0.198 (0.195 assuming Regge behaviour). This value is in good agreement with the 0(a^) 
|3| prediction 0.188 (a s (M|) = 0.109). The non-singlet combination is expected to be less 



sensitive to the low x shape than its singlet counterpart |44[ . Similarly, we observe that the 
value of Ag 3 varies between 1.19 and 1.17 while Au+Ad = Au v + Ad v + 2AS changes from 
0.34 in our model to 0.28 for Regge-type behaviour. We obtain a quite large and positive 
non-strange sea polarization and the whole sea polarization alike (2As = —0.08 seems to 
be reasonable). Finally, AS = 0.26 (0.20) is consistent with existing determinations. 

Performing fits to the inclusive and semi-inclusive data separately we can test the 
impact of each type of data on the total fit. When we use our model to make a fit 
to the data on inclusive spin asymmetries solely, we get x 2 = 95, nearly equal to 96.6, 
which is the contribution of the inclusive data points to the \ 2 = 147 of the total fit 
(x 2 /Nd. .f. = 0.87 is better then in the total case x 2 l^d.o.f. = 0.94). It can be observed 
in Fig.|l], that inclusive spin asymmetries derived from both types of fits do not differ 
much from each other in our model. Moreover results on the integrals over 0.003 < x < 1 
of the structure functions and singlet or non-singlet combinations of quark distributions 
are very close to those obtained in the total fit. In details: AS = 0.27, = 0.126, 

= —0.050, Ag3 = 1.05. Although the total polarization of quarks of a certain flavour 
(i.e. Au(d,s) = Au(d) v + 2Au(d, s)) does not vary much, the division between valence 
and sea quarks differs from the total fit considerably. For example, in 0.003 < x < 1 
region, we get Au v = 0.336, Ad v = —0.565. The reason is, that the valence and sea 
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quark distributions of the same flavour have the same weight (electric charge squared) in 
the inclusive spin asymmetry (eq. (|12"D). Hence the asymmetry is sensitive only to the 
whole Aq distributions but not to valence and sea quark distributions separately. The 
distributions Am and Ad derived from the total fit and the fit to inclusive data have the 
same shape, which can be seen in Fig. |], whereas splits between valence and sea quark 
distributions are different in both cases (compare Fig. |l|). 




Figure 4: The whole Ad and Au distributions derived from the total fit and from fits to 
the data on inclusive and semi-inclusive spin asymmetries separately. 



The comparison to the similar analysis (performed in Ref.0), which uses the MRS 
parametrization as an input shows us the influence of choice of input parametrization 
on distributions and first moments. There is almost no difference between first moments 
of distributions for quarks of a certain flavour obtained in Ref. || and in this analysis. 
Corresponding results are: Am = 0.76 (0.70 assuming Regge behaviour), Ad = —0.52 
(-0.37), As = — 0.07(-0.07). There is a significant difference in division between valence 
and sea quarks. What is most important, the total polarization of the sea quarks changes 
its sign. In Ref.|| Ag sea = — 0.18(-0.22) whereas we have obtained Aq sea = 0.22(0.13). 

Parameterizations obtained using inclusive and semi-inclusive data give a good de- 
scription of the semi-inclusive asymmetries, as can be seen in Fig. |3] (the contribution 
of 48 semi-inclusive data points to \ 2 = 147 of the total fit amounts 50.4, where the 
data point for A\ + at x = 0.48 gives the biggest part). Still, fit performed using only 
semi-inclusive data leads to different set of coefficients of distribution functions and other 
integrated results, as is seen in Figs. |IJ, |4j This is mainly caused by the data point for A\ + 
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at x — 0.48, which makes that semi-inclusive asymmetries go below the total and inclusive 
predictions (Fig. |2|) . Absence of this point would improve the result of comparison to the 
experiment and change high x behaviour of the fitted semi-inclusive asymmetries. 

We get for the fit with semi-inclusive data points only x 2 — 39.3 (x 2 /^d.o.f. — 1-2). 
The integrated quantities are: 

Au v = 0.68 ±0.01 , Am = 0.02 ±0.05 , Am = 0.70 ± 0.05 , 

Ad v = -0.32 ± 0.05 , Ad = -0.02 ±0.07, Ad = -0.34 ± 0.09 , 

As = -0.035 ±0.004 , (22) 

r? = 0.137 ±0.003 , = -0.042 ±0.018 , Tf = 0.044 ±0.010 , 

AS = 0.30 ±0.04 , Ag 3 = 1.08 ±0.21 , Aq sea = -0.06 ± 0.17 . 

There is no important difference between above first moments and integrals obtained 
assuming Regge behaviour for low x. These results are in a good agreement with exper- 
imental estimations |2"B| . The presence of various weights in the semi- inclusive spin 
asymmetries (eq. (|12|)) induces that division of Am and Ad between valence and sea parts 
is no longer strongly model dependent. Also differences among sea quarks of different 
flavours are emphasized. In Fig. ffl, one sees that the parametrizations obtained using 
only semi-inclusive data points give inclusive asymmetries too far from experimental data 
points. If we compute x 2 f° r all data points of both types with the obtained distributions 
we get 505 what is an unacceptable value. The substantial part comes from the El 54 
data for the neutron target, mainly due to the differences in the whole Ad distribution 
obtained from semi-inclusive and total fits, as can be seen in Fig. f|. 

We have performed an analysis of the world data on polarized deep inelastic scattering, 
inclusive and semi-inclusive, assuming that Am ^ Ad, i.e. A5 ^ 0. But in the inclusive 
case only whole quark distributions of a certain flavour are distinguished, i.e. have different 
weights in the asymmetry. As AM(<i) = Au(d) v + AS =F AS, putting AS ^ gives two 
additional coefficients, very weakly constrained. We have performed also a fit to the 
inclusive data only, putting AS = and we have got almost the same value of x 2 as before. 
Although it is possible to obtain the information about difference between Am and Ad 
without taking into account semi-inclusive data, the results are very poorly constrained 
by the data. 
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The situation is better in the semi-inclusive case where all of the distributions (valence 
and sea separately) appear in the semi-inclusive spin asymmetry with different weights. 
Hence, up till now, there are 24 data points^ constraining the coefficients of the AS 
distribution. Performing a fit to the semi-inclusive data with AS = we have got a 
slightly worse value of \ 2 then in the case with AS 7^ 0. 

The inclusion of available semi-inclusive data to the analysis of the inclusive events 
gives more stable results. Using data on semi-inclusive spin asymmetries we can distin- 
guish valence and sea quarks distributions of the same flavour as well as Au and Ad. 
However parametrization obtained as the best fit to the semi-inclusive data gives the 
unacceptable description of the inclusive ones, mainly due to the differences in the Ad 
distribution. Hence, we have got no perfect consistence of inclusive and semi-inclusive re- 
sults in our model. Additional data from semi-inclusive experiments using 3 He target can 
reverse the situation. The next step in the analysis, i.e. addition of the Q 2 -dependence of 
distribution functions can also improve an agreement of our model with an experiment. 

Our analysis shows that the result which gives the polarization of the sea quarks 
depends strongly on used parametrization of the polarized parton distributions. 



2 for the deuteron target AS does not appear in the formula for the semi- inclusive spin asymmetry. 
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